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We study the effect of electrostatic interactions on the distribution function of the end-to-end 
distance of a single polyelectrolyte chain in the rodlike limit. The extent to which the radial 
distribution function of a polyelectrolyte is reproduced by that of a wormlike chain with an adjusted 
persistence length is investigated. Strong evidence is found for a universal scaling formula connecting 
the effective persistence length of a polyelectrolyte with its linear charge density and the Debye 
screening of its self-interaction. An alternative definition of the electrostatic persistence length is 
proposed based on matching of the maximum of the distribution with that of an effective wormlike 
chain, as opposed to the traditional matching of the first or the second moments of the distributions. 
It is shown that this definition provides a more accurate probe of the affinity of the distribution to 
that of the wormlike chains, as compared to the traditional definition. It is also found that the length 
of a polyelectrolyte segment can act as a crucial parameter in determining its elastic properties. 

PACS numbers: 82.35.Rs, 87.15. La, 36.20.-r, 82.35.Lr 



I. INTRODUCTION AND SUMMARY 

The close connection between the elasticity of rodlike 
cytoskeletal polymers and the mechanical properties of 
cells has been extensively documented Q, y, |^ . For ex- 
ample, actin filaments are known to play a crucial role 
in the recovery of eukaryotic cell shape in the face of the 
stresses imposed by cell movement, growth and division. 
Furthermore, the stiffness of DNA constrains and con- 
trols both its storage in the cell nucleus and access to 
it by the proteins and enzymes that are central to its 
role in biological processes. However, there is, as yet, no 
complete, quantitative understanding of the relationship 
between the elastic attributes of these polymers and the 
influences that control them. Such an understanding is 
crucial to the development of a general description of the 
functionality of polymers. 

An additional spur to an increased focus on the con- 
nection between flrst-principle energetics and the elastic 
characteristics of polyelectrolytes is the fact that experi- 
mental techniques have been developed that allow for the 
imaging of single fllaments in solution 4]. This means 
that it is now possible to test models and theoretical 
predictions at the level of a single chain polymer. 

The key challenge in this area is to produce a theoret- 
ical analysis that connects the electrostatic interactions 
and local energetics of a polyelectrolyte with its overall 
mechanical properties. In the case of neutral chain poly- 
mers, the wormlike chain (WLC) model of Kratky and 
Porod 5] provides a powerful and convenient characteri- 
zation of flexibility, through the persistence lengthy which 
is the decay length of the tangent-tangent correlation 
function of a stiff, inextensible polymer chain. The lack 
of a similarly inclusive and readily implemented model 
for polyelectrolytes has been addressed through the ap- 



plication of the WLC model to polyelectrolytes, with the 
introduction of an effective, "electrostatic," persistence 
length 0, • The shortcomings of this notion, which re- 
places the many length scales in a polyelectroyte by a 
single effective persistence length, have been commented 
on H, . Additionally, codifying electrostatic inter- 
actions in terms of single length ignores the possibility 
of chain collapse resulting from fluctuation-induced self- 
attraction 

If a suitably modifled WLC model did, indeed, provide 
a working model for the calculation of the conformational 
properties of polyelectrlytes, then a number of powerful 
results could be brought to bear on the study of this 
system. Among the more useful of these is the recently 
derived end-to-end distribution of short segments of 
an inextensible polymer The distribution, which is 
expected to apply to segments whose backbone length 
is less than the persistence length of the polymer — the 
rodlike limit — was shown to be accurate both in that 
limit and somewhat outside it. The distribution possesses 
the virtue that it can, in principle, be utilized in the 
analysis of currently feasible experiments on polymers in 
the rod-like limit 13]. One can, then, determine whether 
or not the WLC model accurately reflects the properties 
of the polymer segments in question. If the model has 
been shown to apply, then the flt of the distribution to 
the data directly yields the persistence length and, hence, 
the elastic modulus of the polymer. 

An additional advantage of the end-to-end distribution 
function is that it is a function, rather than a number. 
That is, the end-to-end distribution function contains, at 
least in principle, an inflnitely greater amount of infor- 
mation than can be encoded in a single quantity. The 
distribution thus represents a more comprehensive char- 
acterization of the elastic and thermal characteristics of 
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a charged semifiexible chain. In this sense, it provides 
one with an opportunity to disentangle the actions of 
electrostatic and mechanical interactions as they affect 
the conformational properties of a fluctuating polyelec- 
trolyte. 

In this article, we calculate the end-to-end distribution 
of a charged, inextensible, semifiexible chain. This distri- 
bution is then utilized to address two central questions 
relating to the conformational statistics of a polyelec- 
trolyte. Those questions are: 

1. Under what circumstances is the end-to-end dis- 
tribution of the PE reproduced by the end-to-end 
distribution of a WLC with an adjusted persistence 
length? 

2. Under what circumstances will the effective persis- 
tence length of a PE be as predicted by existing 
formulas? 

As will be shown below, we are able to identify 
regimes associated with strong electrostatic coupling, 
weak screening, or high intrinsic flexibility of the PE, in 
which the end-to-end distribution deviates from that of 
a WLC. We are also able to identify the regions in which 
this approximation is accurate. The concept of electro- 
static persistence length is examined, and an alterna- 
tive definition of this important length scale is suggested. 
Our newly-defined electrostatic persistence length is com- 
pared with a well-known and widely-utilized formula for 
the persistence length of a stiff, charged rod 6]. The 
new electrostatic persistence length is found to obey a 
universal crossover scaling law. 

The rest of the paper is organized as follows: Sectionllll 
describes the model that we use to study the elasticity of 
stiff PE's. This model is a straighforward extension of the 
Kratky-Porod model of a WLC. The energy of this sys- 
tem consists of the mechanical energy required to induce 
a local curvature in the chain and the screened electro- 
static repulsion between the charges on the chain, which 
are assumed to be uniformly distributed along it, in the 
form of a constant linear charge density. The distribu- 
tion function of the end-to-end distance of this chain is 
then introduced and given a precise mathematical form. 
An expression for this distribution is derived, under the 
assumption that the chain does not undergo too much 
of a distortion from it minimum-energy configuration as 
a straight line. That expression, which is based on the 
eigenvalues of the Hamiltonian of the PE forms the ba- 
sis of all the results reported in this article. Details of 
the expansion of the electrostatic self-interaction about 
the rodlike configuration are presented in Appendix IXI 
Section IIIII is devoted to a discussion on the structure 
of those energy eigenvalues. The influence of the elec- 
trostatic interaction, which is most pronounced at on 
the lower eigenvalues, is discussed. Also discussed are 
the methods and criteria utilized to guarantee the nu- 
merical reliability of the end-to-end distribution. Section 
IIIIBI brieflv summarizes exact results that have been ob- 
tained for the eigenvalues of the Hamiltonian consisting 



entirely of the unscreened Coulomb interaction. These 
results are relevant to the conformational statistics of a 
chain stiffened entirely by such interactions. A more de- 
tailed description of the way in which the results were 
obtained is relegated to Appendix IbI Section HVl presents 
the results for the distribution function in various situ- 
ations. Particular attention is paid to the question of 
the relationship between the distribution function calcu- 
lated here and the corresponding end-to-end distribution 
function of a neutral rod-like semifiexible polymer. Sec- 
tion scrutinizes the concept of electrostatic persistence 
length. The point of comparison is a general formula 
for the persistence length, introduced by Odijk. We dis- 
cuss the ways in which a persistence length can be ex- 
tracted from the end-to-end distribution, and we justify 
the method of choice in this investigation, which is based 
on the location of the maximum of the end-to-end dis- 
tribution. We find that our results are consistent with 
Odijk's formula when the distribution we calculate can be 
made to collapse on the end-to-end distribution of an un- 
charged WLC. Conversely, our results do not agree with 
the Odijk formula when the distribution is inconsistent 
with that of a WLC with adjusted persistence length. We 
find that the effective persistence length of the end-to-end 
distribution function can be described by a "universal" 
formula, incorporating a scaling form and a new, unex- 
pected power law. This formula holds whether or not 
the effective persistence length is described by Odijk's 
expression. Section IVII focuses on this behavior of the 
effective persistence length. Section IVlII discusses issues 
related to the moments of the distribution, and the use 
of those moments to define a persistence length. Section 
IVIIII concludes the paper. 



II. MODEL ELASTICITY FOR 
POLYELECTROLYTES 

Because of the inextensibility of the polyelectrolytes 
under consideration p^4jl^ , we adopt Kratky and Porod 
wormlike chain model y\ to describe the bending energy 
of the chain. In this model, polymers are represented 
by a space curve r(s) as a function of the arc length 
parameter s. The total energy of the chain, which is the 
sum of the intrinsic elasticity and the electrostatic energy 
can be written as 

knT 2 Jo \ ds J 2 Jo |r(s)-r(s')| 

(1) 

where t is the unit tangent vector. We do not take into 
account the fluctuations in the charges localized to the 
chain and in the counterion system that can give rise to 
attractive interactions leading to chain collapse 10]. The 
above equation contains several different length scales: 
(1) the average separation between neighboring charges 
6, (2) the Debye screening length (3) the Bjerrum 
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length Ib = e^/e/ceT, the quantity e being the dielectric 
constant of the ion- free solvent, (4) the intrinsic persis- 
tence length £po, and (5) the total chain length, L. Two 
of these length scales, namely Ib and 6, always appear 
together in the form of the ratio ^b/^^, which is denoted 
by p in Eq. 

We can thus construct three independent dimension- 
less ratios: /3L (a measure of the strength of the elec- 
trostatic interactions), (a measure of the degree of 
screening), and ipo/L (a measure of the intrinsic flexi- 
bility). A comprehensive exploration of parameter space 
entails the examination of the effects associated with a 
change in the electrostatic coupling and the degree of 
screening. We focus on inextensible chains that are sub- 
stantially stiffened by their bending energy {£po ~ L) and 
for which the bending energy alone is insufficient to keep 
the chain nearly rod-like {£po <C L). The chain is as- 
sumed to be sufficiently stiff that excluded volume does 
not play a role. 

The end-to-end distribution function is defined as fol- 
lows: 



Gir) = {S{r - R)), 



(2) 



where R = r(L) — r(0). The average in Eq. Q is over 
an ensemble of PE chains. The function ^(r) is, then, 
the probability that a given chain in the ensemble will 
have an end-to-end distance equal to r. We make use of 
the procedure that Wilhelm and Frey have implemented 
to calculate the end-to-end distribution function for in- 
extensible neutral polymers pj(. In order to calculate 
^(r) in the vicinity of the rodlike limit we restrict our 
consideration to regions in which the combination of in- 
trinsic stiffness and repulsive strength of the Coulomb 
interaction keeps the chains in their rodlike limit [TEj . 
Using the integral representation of the S function, we 
can write Q{r) as 



X exp • r — z/i • J ds t{s) — J ds ^ 



dt{s) 
ds 



dsds' 



j-«|r(s)-r(s')| 

|r(s)-r(s')| 



(3) 



in which ^ = / Vt{s) exp (— TY/ZcbT). We then parame- 
terize the unit tangent field as 



t(.) 



{aa;{s),ay{s), 1) 
^l+a2(5)+a2(s)' 



(4) 



in which the constraint of inextensibility is manifestly 
satisfied. 

Making use of the relation t{s) — r(s') = du t{u)^ 
we expand the screened Coulomb interaction about 
the rodlike configuration to quadratic order in Si{s) = 
{ax{s), ay{s)) (see Appendix IXI for details). We then use 

the series expansion 3.{s) = a/2 Xl^o (^T^) ^P" 

propriate for the open-end boundary condition, and as- 
sume for simplicity that r is, on average, oriented along 
the z-axis so that Aq = ^ ds a{s) = 0. Substituting 
Eq. (PT|l in Eq. ©, we obtain 



X exp < i/i • r — i/HzL 



2 

n=i 

n=l n,m=l / 



where 



Enm = [ ^{kl^ n^) In 



X < COS 



COS 



-{n — m) 



TT 



(n + m) 



sin (^) sin (^-^ 

kz [kz - (n - m)7i/L\ 

"sin (i^) ,m (ti - tol) 



kz[kz - (n + m)7r/L] 
(7r/6)2 + /t2 + fc2- 



sin - 2f ) sin ( 



kzL Trm\ 

2 2j 



(kz — nn / L){kz — mn / L) 



sin + ^) sin ( 



^in^ (¥) 

k^ 



r 



kzL Trm ^ 

2 2) 



{kz + nn / L)(kz — mn / L) 



(6) 



are the elements of the electrostatic energy matrix in the cosine basis set, with 6, the average separation between 
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neighboring charges, serving as the short distance cutoff. 

Performing the functional integral over A^,, we then 
obtain 



C?(r) = 



2^ 



,ia;(l-r^/L) _ 



det [T] 



det [T + iujl] 

ij) 

where T = {rni)'^ / L)l + jSLE, I is the identity ma- 
trix, and uo — —fizL. Note that the distribution func- 
tion in Eq. Q is manifestly normalized to unity, i.e. 
/ d^r Q{r) = 1. We can then integrate over the orient a- 
tional degrees of freedom to find the radial distribution 
function 



G(r) 



(8) 




n 



FIG. 1: Eigenvalues of electrostatic energy matrix calculated 
from Eq. at L/b = 1000 and kL = 10 (dashed curve), 
compared to the eigenvalues calculated in the infinite length 
limit in Eq. (|14|) (solid curve). 



that is subject to the normalization dr G{r) = 1. 
Using r = (r sin cos (/), r sin sin (/), r cos 0) , we obtain 



^iujil-r/L) _^ giu;(l+r/L) 



det [T] 



det [T + iul] 
(9) 



which can be rewritten as 



G(r) 





duo 




J — OO 






oo > 

n 

n=l 




.An 4 


- iuj J ^ 



,iu;(l-r/L) ^ gia-(l+r/L) 



(10) 



in terms of the eigenvalues. An, of the matrix T. Contour 
integration then yields 



-A„(l-r 



with 



fin) = n 



A,; 



An ~ Aj 



/L) 
(11) 

(12) 



In the above expressions, the second term in brack- 
ets contributes negligibly in the rodlike limit, because 
even the smallest eigenvalues are considerably larger than 
unity in this regime. Neglecting this term thus yields 



G(r) = mL -r)J2 ^nf{n) e-^-^i"'' 



/L) 



(13) 



where the normalization prefactor JV is used to compen- 
sate for the small error caused by neglecting the second 
contribution. This prefactor can be found by enforcing 
normalization on G(r). 



III. ENERGY EIGENVALUES 
A. General Results 

To obtain the radial distribution function we must de- 
termine the spectrum of the effective Hamiltonian T. 
Since this effective Hamiltonian incorporates the energy 
cost when the PE deviates from a rodlike configuration, 
the spectrum of its eigenvalues An should be positive and 
a monotonically increasing function of the index n. Note 
that it is the lowest eigenvalues that exert the dominant 
infiuence on the radial distribution function. 

Figure ^ shows the eigenenergies of the electrostatic 
interaction E at i^L = 10. To emphasize the importance 
of the finite length L of the PE, we calculate the matrix 
Enm [using Eq. Q] in the limit L ^ oo, and find 



Enm{L oo) 



In 



(n7r)2 + (kL)2 



(n7r)2 
(utt)^ + (kL)^ 



1 . (14) 



This limit can be achieved in Eq. (|T^ by us- 
ing limL^oo sinjL (^ - ff )] / - ff ) = 2T:5{k, - 
nir/ L). Figure Qalso displays the eigenvalues of the diag- 
onal energy matrix in Eq. (|T^ . A comparison of the two 
curves clearly highlights the importance of end effects for 
PE chains of finite length. 

The electrostatic energy comprises one portion of the 
total energy of the semifiexible PE. The other contribu- 
tion arises from the elastic modulus, quantified by the 
"bare" persistence length, ipQ and appearing in the en- 
ergy as the first term on the right-hand side of Eq. (QJ. 
The large eigenvalues of the matrix T are dominated 
by the diagonal terms {nn)'^ {£po / L) in, for instance, the 
right-hand side of Eq. (|5J. Thus, the effect of electro- 
static interaction is more pronounced on the lower eigen- 
values, while it is swamped by semifiexible energetics at 
short length scales. In Fig. [21 the eigenvalues of the 
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of matrices by looking for convergence of the numerical 
procedure, as noted above. We then check to see that the 
size, of the basis set satisfies N < L/b. At no point 
in our calculations was this inequality violated. 

As r ^ 1, the dimension of the matrices must be very 
high in order to obtain reliable answers. It is possible 
to calculate the integral in Eq. (|Tn|) numerically with- 
out performing contour integration. This method is very 
time-consuming but it is more reliable because this way 
we reduce significantly the roundoff error caused by the 
sum in Eq. (|T^ . We used this approach for very large 
values of PL. 



FIG. 2: Eigenvalues of an elastic charged rod with ipo/L = 
0.5, kL = 0, and pL = 500. The bending eigenenergies are 
distinguishable from total ones only for lower values of n. 
Comparison of the first 20 eigenvalues of bending (solid line) 
and total (dashed line) eigenvalues could be seen in the inset. 



total Hamiltonian are compared with the bending eigen- 
values. Only the first few eigenvalues are distinguishable 
from each other. The value of n for which the electro- 
static effects can be neglected depends on the strength of 
electrostatic interaction PL and the intrinsic persistent 
length. 

The general behavior of the larger eigenvalues associ- 
ated with "pure" electrostatic energy, which grow log- 
arithmically with index, n, is fundamentally different 
from those resulting to "pure" elastic energy, which grow 
quadratically with index. In the case of a neutral poly- 
mer with an intrinsic persistence length that is compa- 
rable with the total length of polymer, the series in Eq. 
(|T^ converges very rapidly, as noted in Ref. . When 
the chain is charged, and especially when the persistence 
length is due predominantly to electrostatic effects, more 
terms in the series Eq. ([T^ must be preserved in order 
to obtain a stable answer for G{r). In our calculation, 
the matrix T is truncated at a size much greater than 
the number of terms needed to obtain an accurate an- 
swer for the series in Eq. (|13|) . There are two reasons for 
this. First, the greater the dimension of the truncated 
matrix, the more accurate are the values of lower eigen- 
values which participate in the sum. Second, retaining 
more terms in the product on the right-hand side of Eq. 
(|T^ helps us achieve a higher accuracy in our result for 
the function f{n). In practice, we increased the dimen- 
sion of the matrices until the effect (on the end-to-end 
distribution function) of a further increase in the size of 
T by a factor of four was less than a part in a thousand. 

The truncation of the matrices above is also related 
to the short distance cutoff b. The dimensions of these 
matrices should not be larger than L/b. This is because 
matrices corresponding to a cosine basis set for which 
the wavenumber is greater than l/b are infiuenced by 
features at length scales shorter than 6, and such fea- 
tures are inconsistent with the inherent coarse-graining 
in our energy expression. We determine the truncation 



B. Exact Analytical Eigenfunctions and 
Eigenvalues of the Unscreened Coulomb Interaction 
Energy 

Although we are not able to produce exact, analyti- 
cal solutions to the eigenvalue equations of the quadratic 
weight in Eq. (0), for arbitrary values of the inverse 
screening length, it has proven possible to obtain the 
eigenfunctions and eigenvalues of the unscreened version 
of the energy operator. The results of the analytical in- 
vestigation provide a check for our numerics and provide 
insight into the qualitative behavior of eigenfunctions and 
eigenvalues of the energy operators controlling the con- 
formational statistics of rod-like PE's. 

As it turns out, the investigations of the properties of 
the operator in this special case is most conveniently car- 
ried out in "real space" — that is, with arc- length, 5, as 
the independent variable. In this section, we will summa- 
rize the methods utilized and the results obtained in this 
investigation. The details of the calculation are presented 
in Appendix El The key result is that the eigenfunctions 
of the unscreened Hamiltonian have the form 

=s-'/'sin(fcln^+0(fc)) (15) 

and that the eigenvalues are given by 

A(fc) = A(^-i + ifc^ (16) 

where the function A{p) is defined in Eq. (|B12|) . When 
the dimensionless wave- vector-like quantity, /c, is large, 
the eigenvalues are essentially linear in \nk. This is to 
be expected, given the form of the unscreened Coulomb 
interaction, and is, in fact, consistent with our numerical 
results. 

We can now attempt a comparison between the expres- 
sions for the eigenfunctions in Eq. (|T5j) . and eigenvalues 
in Eq. (|T7)j) . and the numerically generated eigenfunc- 
tions and eigenvalues of the unscreened energy operator 
whose matrix elements in the cosine basis are given in 
Eq. 0. The eigenvalues are compared in Fig. 01 and a 
comparison for a low-lying eigenvector is displayed in Fig. 
^ As indicated by the figures, the numerical correspon- 
dence is excellent. The quality of the comparison serves 
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as a validation of the numerical calculations carried out 
here. 




index 

FIG. 3: Comparison of the eigenvalues A(A:) as given by 
Eq. (|B14|) with the numerical eigenvalues of the unscreened 
Coulomb energy matrix as given by Eq. 0. As indicated 
in the figure, the horizontal axis is the index of the eigen- 
value. The formula (fB14|) is the solid curve, passing through 
the filled circles, which indicate the location of the eigenvalues 
of the discrete version of the operator in Eq. To facil- 

itate the visual comparison, the numerical eigenvalues have 
been "thinned out," in that every third numerical eigenvalue 
is shown. 



resulting from an unscreened Coulomb interaction 

^bending ^ IpO [(1 + 4fc^)/4] [fcV(l + fc^)] 
^Coulomb 2lB\{k)HL/b) ^ ' 

Note that in the limit of an infinitely long polyelectrolyte 
chain, the unscreened Coulomb interaction always dom- 
inates. The divergence in the denominator as L ^ oo 
is, however, quite slow. As practical matter, the bending 
energy due to the intrinsic cost of introducing curvature 
into a semi-flexible chain will eventually overcome the 
Coulomb energy as a contribution to the conformation 
energy of the chain. In light of this fact, one can obtain 
results for the end-to-end distribution of an unscreened 
chain in the presence of a bending energy. 

In this regard, it should be noted that in the absence 
of a bending energy, the expression (|10|) for the end-to- 
end distribution of an unscreened polyelectrolyte chain 
is formally divergent. This is because of the very slow 
increase in the eigenvalues ([T^ as a function of k. A 
full investigation of the conformational statistics of the 
unscreened polyelectrolyte chain has not yet been carried 
out, and will in all probability yield surprises. 

IV. RADIAL DISTRIBUTION FUNCTION OF 
POLYELECTROLYTES 




-5 



-10 
-15 

FIG. 4: Comparison of the numerically determined eigenvec- 
tor of the discrete version of the Coulomb energy kernel asso- 
ciated with the sixth-lowest eigenvalue (solid curve) and of the 
eigenfunction generated by the formula (|15f) (dashed curve). 
The value of the parameter k in that formula is adjusted to 
optimize the fit between the two curves, with the use of a 
least-squares procedure. The actual value of k used here is 
1.3. The oscillations in the solid curve are the result of the 
cutoff in the cosine basis set. 

It is possible to assess the importance of a small term 
representing the influence of "pure" bending energy, by 
treating that term as a perturbation of the overall energy 
of the polyelectrolyte chain. At first order in perturba- 
tion theory, a straightforward calculation leads to the 
following result for the ratio of the pure bending energy 
contribution to the eigenvalue to the zeroth order energy 



Using the eigenvalues of the effective Hamiltonian for 
the rodlike PE, we can calculate the radial distribution 
function from Eq. ([T^ . Figure [51 illustrates the effect 
of the screened electrostatic interaction on the distribu- 
tion function of a neutral chain (dashed line), for three 
different values oi nL = 12, nL = 6 and a^L = 3. 
/3L = 100 for the three PE distributions. As mani- 
fested in the figure, the distribution function of the poly- 
electrolyte has shifted toward larger extensions and is 
peaked more sharply around the maximum compared to 
the neutral chain. Upon decreasing the screening length, 
the monomer-monomer repulsion becomes more strongly 
screened and the distribution function of the end-to-end 
distance peaks towards shorter extensions. 

To construct a quantitative measure for the elasticity 
of polyelectrolytes, we investigate to what degree and 
under what conditions the radial distribution function of 
a polyelectrolyte can be reproduced by the distribution of 
a neutral chain with an adjusted persistence length. To 
this end, we attempt to collapse the distribution function 
of a polyelectrolyte onto the distribution function of a 
wormlike chain with an effective ip^ and observe that the 
answer to the above question depends on the values of 
ipo/L^ kL and [3L. 

We investigate the effects of overall charging of the 
chain as quantified by the combination of the screen- 
ing length, equal to of the length of the chain, L 
and of the chain's intrinsic stiffness, described in terms 
of the bare persistence length, We are able to iden- 
tify regimes in which the conformational statistics of a 
rod-like PE is effectively indistinguishable from that of 
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FIG. 5: Comparison of distribution function of polyelec- 
trolytes for a^L = 12, a^L = 6 and kL = 3, respectively from 
left to right. 

a neutral WLC, and we also find that there are regions 
in parameter space in which the radial distribution of a 
charged semiflexible rod-like chain cannot be reproduced 
by the corresponding distribution of a neutral chain. 

Our conclusions are summarized in the four subsec- 
tions below. Briefly, we find that sufficiently strong 
charging, a sufficiently long screening length, sufficiently 
short overall chain length or sufficiently weak intrinsic 
stiffness can give rise to deviations between the confor- 
mational statistics of a PE and that of a WLC. 



A. Effect of Electrostatic Charging 

There is a substantial regime in the parameter space 
in which the distribution of a polyelectrolyte matches 
exactly that of a wormlike chain with an adjusted per- 
sistence length, particularly when the electrostatic inter- 
action plays a perturbative role in the chain energetics. 
There are also regimes in which the electrostatic inter- 
action gives rise to a substantial portion of the stiffening 
energy, and in which the radial distribution of the PE is 
the same as that of a neutral WLC. For a sufficiently 
strong electrostatic interaction, however, the polyelec- 
trolyte distribution begins to deviate from the wormlike 
chain form. This is illustrated in Fig. [SI where the poly- 
electrolyte distribution is compared with that of the best 
effective wormlike chain description for the two cases of 
low and high charging. It is clear that at strong enough 
charging, the electrostatic energy establishes the confor- 
mational statistics. 



B. Effect of Electrostatic Screening 

The degree of screening of the electrostatic interac- 
tion similarly affects the form of the end-to-end distribu- 
tion. In the high-ionic-strength regime it is possible to 
obtain a satisfactory match of the polyelectrolyte distri- 
bution with that of an effective wormlike chain. How- 
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FIG. 6: (a) The distribution for a polyelectrolyte with ipo = 
0.5, PL = 600, and kL = 10, compared to that of a wormlike 
chain with ip = 1.3. Inset: (b) The distribution for a poly- 
electrolyte with ipo =0.5, (3L = 60, and kL = 10, compared 
to that of a wormlike chain with ip = 0.6 (two superimposed 
curves) . 
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FIG. 7: The distribution for a polyelectrolyte with IpQ = 0.5, 
pL = 600, and kL = 50, compared to that of a wormlike 
chain with Ip — 0.56 (two superimposed curves). 



ever, as screening decreases, the polyelectrolyte distribu- 
tion deviates significantly from the wormlike chain form. 
Comparison of Figs. El^a) and illustrates the differ- 
ence between high and low screening for a given intrinsic 
stiffness and charging. We observe that as we increase 
screening, the two distribution collapse on top of each 
other. Our calculations verify that under physiological 
conditions = 1 nm~^), the distribution functions for 
rod-like DNA segments (L^lOO nm) as well as those of 
stiffer actin filaments also collapse onto the end-to-end 
distribution for neutral chains with an effective persis- 
tence length given by Eq. ([T^ . 



C. Effect of Finite Length 

In most of the studies of the elasticity of polyelec- 
trolytes that have been carried out to this point, the finite 
length of the polyelectrolyte and the corresponding end- 
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FIG. 8: (a) The distribution for a polyelectrolyte with 
^po/L = 0.0001, pL = 36000, and kL = 100, compared to 
that of a wormUke chain with ip = 0.87 (two superimposed 
curves). Inset: (b) The distribution for a polyelectrolyte with 
ipo/L = 0.01, f3L — 360 and — 1, compared to that of a 
wormlike chain with — 0.565. The difference between the 
two sets of parameters is only in the length of the chain. 

effects have been assumed to be unimportant. We have 
re-examined this assumption by considering two sets of 
parameters that correspond to identical values for -^^o, 
/3, and and different values for the contour length L 
of the polyelectrolytes. As can be seen in Fig. |S1 while 
the distribution for longer chains can be satisfactorily 
collapsed onto that of a wormlike chain, the situation is 
completely different for shorter chains. This highlights 
the importance of end effects in the elasticity of polyelec- 
trolytes. 

It is important to note that a rescaling of the backbone 
length of the PE is not an acceptable stratagem for im- 
proving the agreement between the PE radial distribution 
and that of a WLC. This is because the backbone length 
is essentially fixed by the rod-like chain condition. The 
shortening and thickening that is associated with inter- 
mediate blob-like structures will not occur. In fact, 
we have observed that a reduction in the effective value 
of L actually degrades the quality of the correspondence 
between the conformational statistics of a PE and those 
of the corresponding WLC. In general, our calculations 
indicate that the difference between the radial distribu- 
tion of the PE and the WLC is partly due to the influence 
of end effects. A complete investigation of these effects 
is described in Ref. [TJ- 



D. Effect of Intrinsic Stiffness 

Fig. IHIJb) indicates that for nL = 1 and jSL = 360 
when the intrinsic persistence length of a chain is small 
{ipQ <C L = 0.01), the radial distribution function of 
a polyelectrolyte does not exhibit near-perfect matching 
with that of an uncharged wormlike chain . In this al- 
most unscreened case, matching is achieved, when the 
chain becomes sufficiently stiff. If ipo/L is increased to 
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FIG. 9: The diagram delineating the regimes in which the 
statistics of a rod-like PE is the same as that of a WLC with 
a suitably modified persistence length from the regimes in 
which the statistics of the two models are fundamentally dif- 
ferent. The curves separating the two regimes correspond to 
^eZ-^Odiik = 0.58 for various values of ipo/L. These lines also 
indicate the values of poL in Eq. (fll^ at different values of 
kL and ipo/L. No curves were drawn outside of the range of 
rodlike behavior. For this reason, the curves for ipo/L = 0.1 
and ipo/L — 0.3 are truncated. 

35, while tzL remains equal to one and /3L = 360, the 
radial distribution of a PE will be indistinguishable from 
the corresponding distribution of a neutral WLC with 
a persistence length ip = 39 L. The fact that such a 
large intrinsic persistence length is required to achieve 
this kind of matching of the two distribution highlights 
the importance of electrostatic interactions in this case. 



E. Applicability of the WLC Model: A General 
Diagram 

The effects discussed above can be summarized in 
a general diagram in the three dimensional parameter 
space. This diagram, shown in Fig. El illustrates the 
conditions under which a WLC model provides an ac- 
curate description of the conformational statistics of a 
short, rod-like PE. The lines shown indicate the loca- 
tions of the crossover regions that separate the regime 
in which one can think of a short segment of a PE as a 
WLC with a modified persistence length from the regime 
in which the statistics of a PE is fundamentally distinct 
from that of a WLC. 

As can be seen in Fig. |H1 for a given value of the screen- 
ing parameter, charging up the PE will cause a crossover 
from WLC into non-WLC behavior. This crossover is, 
however, hampered by increased intrinsic stiffness. 



V. ELECTROSTATIC PERSISTENCE LENGTH 

In this Section, we investigate the concept of the elec- 
trostatic persistence length. As in the previous section. 
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we will find that there are regimes in which Odijk's for- 
mula [Eg. (|T^ below) for the effective persistence length 
of a PE Q is accurate, and that there are regimes in 
which it does not hold. As it turns out, there is a strong 
correlation between the accuracy of this formula and the 
correspondence between PE and neutral WLC conforma- 
tional statistics. When the radial distribution function 
of a PE is well-reproduced by that of a neutral WLC, 
the persistence length of the corresponding WLC is ac- 
curately predicted by Odijk's formula. Conversely, when 
the two distributions do not collapse on one another, the 
persistence length of the "best-fit" WLC is not predicted 
by that formula. 

We have seen in Sec. IIVI that one can collapse the 
distribution function of a polyelectrolyte onto the dis- 
tribution function of a wormlike chain with an adjusted 
persistence length whenever the Coulomb interaction is 
no more than a perturbation to the mechanical stiffness of 
the chain, or when Debye screening is sufficiently strong. 
In these regimes, the persistence length of the neutral 
chain follows Odijk's prediction, in that, = + ^po, 
where is the effective persistence length of the charged 
chain, and the electrostatic persistence length ^e is given 
as ^] 
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(18) 



which reduces to ^osf = P/^ti^ for large /^L 1^0. For 
instance, at [3L = 600 and hzL = 50, if we add ^osf = 
0.06 to ipo = 0.5, we find ip = 0.56. The distribution 
function of this polyelectrolyte collapses perfectly onto 
that of a wormlike chain with = 0.56. It is important 
to note that the expression for ^odijk was derived under 
the assumption that £po ~ the contour length of the 
chain, L, is of the order of its intrinsic persistence length, 
£pO and that the electrostatic interaction has the limited 
effect of "perturbing" the wormlike chain shape of the 
chain 6]. Our results confirm the validity of OSF formula 
in the regime in which it is expected to be correct {hcL ^ 

For long chains, substantial charging is required in or- 
der to enforce the rodlike limit. It is important to note 
that despite the substantial charging of the chain, the 
ratio of the length of the PE to the screening length, /^L, 
must be sufficiently large. This requirement is essential in 
order to minimize the strong influence of the end effects 
on conformational statistics of charged chains. 

The important characteristic of this regime is that 
£pO ^ £e and thus electrostatic energy no longer plays 
a perturbative role. However, we find that in this regime 
Odijk's formula works perfectly well as well, provided 
the screening is strong enough. For example, at £po/L = 
0.0001, pL = 36000, and = 100, we find i^/L = 0.87 
which is a near-perfect match with ^odijk (see Fig. |H|. 




FIG. 10: The intrinsic persistence length of the charged chain 
is ipo/L = 0.5, pL = 45 and k,L = 2. The persistence length 
of the wormlike chain is £p/L — 0.767. 



The accuracy of Odijk's formula when ^ ipo is not 
at all obvious, as OSF was derived in the regime where 
electrostatic interaction plays a perturbative role. 

We emphasize that the reason for the high quality of 
the match with OSF in this regime is different from the 
reason for the corresponding result obtained by Khokhlov 
and Khachaturian 18] for weakly charged flexible chains. 
In our case, the chain is stiff in all length scales; the 
possibility of renormalizing the length and/or charge is 
thus excluded in our formulation. 

As noted in Sec. IIVI increasing the electrostatic cou- 
pling or screening length, one encounters a regime in 
which it is no longer possible to obtain near-perfect col- 
lapse of the polyelectrolyte distribution onto that of an 
effective wormlike chain. FigurenS displays the polyelec- 
trolyte end-to-end distribution (solid curve) along with 
the modified wormlike chain distribution (dashed curve) 
in a case in which it is possible to obtain a good, but not 
perfect fit. The fit was obtained by matching the location 
of the maxima of the two distributions, and the electro- 
static persistence length attributed to the polyelectrolyte 
distribution is that of the wormlike chain associated with 
the dashed curve. The PE distribution shown in the fig- 
ure matches the result of computer simulation of a rod- 
like PE |19|. In this case the ratio of electrostatic persis- 
tence length extracted from the distribution function to 
Odijk's persistence length is, ^e/^Odijk = 0.58. As noted 
at the beginning of this section, the more the two dis- 
tributions are different from each other, the more the ef- 
fective persistence length of the polyelectrolyte deviates 
from the Odijk's prediction. For example. The persis- 
tence length of the wormlike chain shown in Figure [SI (b) 
is £p = 0.565, while the effective persistence length of the 
polyelectrolyte according to Odijk is £po + ^odijk = 4.52. 
The two distributions differ significantly and the ratio is 
4/^Odijk = 0.12. 

Our general observation is that the ratio of £e extracted 
from the effective wormlike chain distribution to Odijk's 
persistence length correlates with the quality of the fit 
of an effective wormlike chain end-to-end distribution to 
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FIG. 11: Comparison of our results for the electrostatic per- 
sistence length with Odijk's finite size formula (solid line) at 
kL = 2 with ipo/L = 0.5 (filled circles) and ipo/L = 0.01 
(hollow circles). The inset is for kL = 17 and ipo/L = 0.5. 



FIG. 12: Comparison of our data, suitably rescaled, to the 
expression in Eq. (|19f) for the form of the correction to Odijk's 
formula for the electrostatic persistence length. The points 
correspond to our rescaled plots of ie/L versus f3L at kL — 
6,8,10,12. 



that of a polyelectrolyte. When this ratio is equal to 
one, the polyelectrolyte can be completely described in 
terms of a wormlike chain. As this ratio decreases, the 
deviation becomes more pronounced. Figure El displays 
a diagram, which delineates the quality of the fit. The 
line that is used in the Figure to separate the two regimes 
corresponds to ^e/^Odijk = 0.58. As indicated in this Fig- 
ure (and described in Sec. II V E() . for fixed nL^ when jSL 
is below a certain value the polyelectrolyte behaves like 
a wormlike chain, while for larger jSL the two distribu- 
tions start to differ significantly. This crossover scale is 
sensitive to the intrinsic flexibility of the PE, as shown 
in Fig. El 

There is a substantial difference between the elec- 
trostatic persistence length of a rodlike polyelectrolyte 
and Odijk's prediction as we increase the electrostatics 
strength I3L. As illustrated in Fig. ^2 the deviation 
of ^e/L from ^odijk/^, with increasing jSL is more pro- 
nounced at lower values of nL. At i<iL=ll^ this deviation 
becomes evident when jSL ~ 1000. However, at kL = 2, 
the deviation is evident already when /3L ~ 30. 

A review of the current literature on elasticity of poly- 
electrolytes reveals that there is no simple theory for 
computing ie with arbitrary intrinsic stiffness. It is gen- 
erally believed that as long as a polyelectrolyte chain is 
stiff, the dependence of £e on tz is correctly given by the 
OSF formula. However, our results indicate that in cer- 
tain regimes the electrostatic persistence length, i^^ de- 
pends on £po, and is smaller than is predicted by the OSF 
formula. Figure ITTl shows that the electrostatic persis- 
tence length of a rodlike polyelectrolyte, ie/L^ depends 
on ipo/L. This dependence is not present in Odijk's for- 
mula. 

Our general observation of the curves of ie versus 
PL is that they asymptote to a power law of the form 
£e (X {PLy where the exponent x{hzL) is less than 1. 
This behavior has been predicted by Barrat and Joanny 
with an exponent x = 1/2 in the limit /3L oo 8]. We 
find, however, that the asymptotic limit is approached 



much more slowly than it is predicted by the Barrat- 
Joanny crossover formula. The section below describes 
our attempt to systematize the relationship of the effec- 
tive persistence length to Odijk's prediction ^201] . 



VI. UNIVERSAL BEHAVIOR 

In search of a possible universal pattern which might 
relate ipo/L^ PL and tzL to ^e, we have rescaled our 
graphs for different values of the parameters. To this 
end, we choose a value of P called Po such that ^e/^Odijk = 
0.58, and rescale the parameter p with it, where the value 
of 0.58 is completely arbitrary. We have followed this 
procedure to rescale our graphs for different values of 
tvL^ keeping ipo/L constant. Figure IT^ contains these 
rescaled graphs for f<iL = 6,8, 10, 12, which collapse sat- 
isfactorily on top of each other. It is important to note 
that when we choose any other ratio of ^e/^Odijk to rescale 
our graphs, we obtain the same universal behavior among 
our curves. The solid line in Fig. IT^ corresponds to the 
following crossover or interpolation formula for ie 



^Odii 



ijk 



■c{p/por' 



(19) 



plotted with the exponent x = 0.4. In the above equa- 
tion, c = 0.724 is constant which ensures the ratio of 
^e/^Odijk for a givcu value of Pq. 

As exemplified in Fig. El this formula provides a re- 
markable fit to our data. To calculate the exponent x, we 
have attempted to fit our curves of log(£odijk/^e — 1) ver- 
sus \og{P/po) to a straight line. As shown in Fig. ^1 all 
our data are fitted to a straight line with slope approxi- 
mately equal to 0.4. Similar curves have been generated 
for other, smaller, values of ipo/L. The behavior of these 
other curves is the same as the one plotted in Fig. 1121 



11 




FIG. 13: Demonstration of the universal scaling behavior for 
le- The points correspond to our rescaled plots oi le/L versus 
f3 L a,t hvL — 6, 8, 10, 12. The equation of straight line is y — 
0.006 + 0.401X. 



The quantity /3o is an increasing function of and 
n. We have been able to obtain the value of jS^ from 
our data. The crossover discussed in the previous section 
occurs when P = Po which corresponds to £e — 0.58£odijk 
in Eq. (^3- Hence, the curves in Fig. El indicate the 
values of /3oL for different values of f<iL and ipo/L. The 
dependence of Po on the parameters £po^ and L seems 
to be in general quite complicated. For Ip^n > 1, we seem 
to find that Pq ~ L'^n^ipo g{nL)^ where ^ is a polynomial 
function. For Ip^n < 1, however, the leading dependence 
on K seems to be a much stronger power law. We are still 
investigating the dependence of Po on the parameters in 
all the different regimes |21|]. 



VII. FIRST AND SECOND MOMENTS OF THE 
DISTRIBUTION 

Much theoretical work on electrostatic persistence 
length, is based on the calculation of (r) or (r^) of a 
polyelectrolyte. Odijk uses the second moment to test 
the validity of OSF when ie > ipo 0. It is thus crucial 
to learn to what extent the conformational statistics of 
polyelectrolytes are similar to those of wormlike chains 
when they share the same (r) or (r^). It is not obvious 
that the matching of first or second moments of the end- 
to-end distribution will yield the matching of the higher 
moments of the distribution as well. 

We extract an "effective" persistence length by match- 
ing the first or second moments of the end-to-end distance 
of a polyelectrolyte distribution with that of an effective 
wormlike chain in order to compare our results with other 
existing theories. 

FigureElcontains the end-to-end distribution function 
for a screened polyelectrolyte with hcL = 17, ipo/L = 0.5 
and pL = 1000 along with the distribution function of an 
uncharged wormlike chain with £p/L = 1.22. The shorter 
distribution is adjusted so that its second moment is the 
same as that of the polyelectrolyte distribution. It is 
noteworthy that the persistence length of the uncharged 
wormlike chain agrees perfectly well with Odijk's predic- 



G(r) 




FIG. 14: The taller distribution (solid line) belongs to a poly- 
electrolyte at kL = 17, ipo = 0.5 and PL = 1000. The shorter 
distribution (dashed line) is adjusted so that its second mo- 
ment is the same as that of the polyelectrolyte distribution. 
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FIG. 15: The thin solid line is the plot of G{r) for a polyelec- 
trolyte with pL = 600, kL = 0, and ipo/L = 0.5. The thick 
solid line with ipo/L = 3.1 belongs to a WLC distribution for 
which its maximum is matched to that of the PE distribution. 
The dashed plot is for a neutral chain with £po/L — 4.86. The 
dashed distribution and polyelectrolyte have the same first 
moment. 



tion, while the two distribution functions are distinguish- 
able from each other. In general, we observe that the 
maximum of the wormlike chain distribution is shifted 
toward larger extension when we match (r^) of the two 
distributions as exemplified in the Figure. 

Similar results are observed when we match the first 
moment of distribution of a polyelectrolyte with that of 
a wormlike chain with an adjusted persistence length. In 
Fig. ^1 the distribution function of a polyelectrolyte 
along with two distribution functions for effective worm- 
like chains is presented. The dashed line and polyelec- 
trolyte share the same first moment. As illustrated in 
the Figure, the two distributions differ significantly. 

We have compared the electrostatic persistence length 
given in Eq. (|T^ with our results obtained by matching 
(r^) of a PE distribution to that of a WLC distribution. 
Figure El shows that for tzL = 17, our data begins to 
deviate from Odijk's finite size formula when PL ~ 1200. 
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FIG. 16: Comparison of our results for the electrostatic per- 
sistence length with Odjik's finite size formula [ Ea. (fTS|) ] (solid 
line) at kL = 17. Our numerical results are based on match- 
ing of the second moment of distributions. 

It is important to note that when we match the max- 
imum of distributions to extract an effective persistence 
length, the deviation from Odijk's prediction begins at 
pL - 400 for = 17 (the inset in Fig. This is 

the point when the two distributions cease to closely re- 
semble each other. However, if we use the matching of 
the first or second moment to calculate the effective per- 
sistence length, the two distributions would clearly be 
distinguishable from each other well within the regime 
in which Odijk's formula for the persistence length re- 
mained accurate. 

This points to the fact that replacing a PE chain by a 
WLC, when they share the same (r^) or (r), is not well- 
justified in all regimes and that one should use care in 
the utilization of the notion of an electrostatic persistence 
length based on matching of first or second moments. 



VIII. CONCLUSION 

The difficulty in producing a complete characteriza- 
tion of the mechanical properties of a charged, semiflexi- 
ble chain arises from the existence of a number of length 
scales in that system. In particular, it is far from clear 
that the notion of a persistence length provides an ade- 
quate description of the mechanical and thermal charac- 
teristics of such a chain. In light of this, the end-to-end 
distribution function provides an attractive alternative. 
Given that the distribution is a function, rather than a 
single number, it represents considerably more informa- 
tion about the charged chain. 

One important application of this distribution is in 
the assessment of the utility of the notion of a persis- 
tence length, in that it is possible to compare the dis- 
tribution obtained experimentally, via simulations, or as 
the result of explicit calculations, with the correspond- 
ing distribution of a neutral semiflexible chain. On the 
basis of explicit calculations, we have been able to deter- 
mine the extent to which the end-to-end charged chain 



can be collapsed onto the end-to-end distribution of a 
neutral worm-like chain with an adjusted persistence 
length. Among the regimes in which this correspondence 
is achieved are those in which electrostatic effects play a 
perturbative role. In addition to collapse of the two dis- 
tributions in these regimes, we also find that the electro- 
static persistence length is given by the formulas of Odijk 
and OSF. This result is consistent with the assumptions 
underlying those formulas. We also observe collapse of 
the distributions and can verify the validity of the Odijk 
and OSF results for the effective persistence length in 
regimes in which the charging of the chain is sufficiently 
strong that electrostatic effects dominate purely mechan- 
ical energetics as long as Debye screening of electrostatic 
interactions is sufficiently strong. 

In fact, we find that there is a strong correlation be- 
tween correspondence of the distributions of charged and 
uncharged chains and the validity of Odijk's formula for 
the effective persistence length. When the distributions 
can be collapsed onto each other, the formula proves to be 
accurate, while lack of correspondence of the two distri- 
butions is accompanied by inaccuracy of the Odijk result 
for the effective persistence length. Our results indicate 
that the difference between the radial distribution of the 
PE and the WLC can be attributed, at least in part, to 
the influence of end effects. In fact, we believe that the 
behavior of the persistence length is substantially con- 
trolled by end effects. One way of understanding this is 
in terms of Odijk's derivation of the expression ([T^ for 
the electrostatic persistence length 6]. This derivation 
is based on a calculation of the energy of a bent segment 
of a charged rod. A key assumption in this derivation is 
that the segment takes the form of an arc of a circle. End 
effects are readily associated with the difference between 
the shape of a real bent rod and the circular arc assumed 
in Odijk's derivation. An exploration of these effects in 
this context is described in Ref 17]. 

Another important finding is that an effective persis- 
tence length, obtained by locating the maximum of the 
distribution, can be described in terms of a scaling for- 
mula, Eq. (|T^ . This formula relates the actual persis- 
tence length to the Odijk predictions. The formula is 
"universal," in that it has a general form that is inde- 
pendent of the parameters utilized, and it incorporates a 
power law that does not appear to be anticipated in the 
Hamiltonian governing the system, nor does it arise from 
any simple dimensional analysis. At this point, we have 
no explanation for either the universal form or the power 
law. 

As noted above, the effect of counter ion condensation 
has been ignored throughout the above work. It has been 
shown that counterion condensation modifies the bend- 
ing rigidity of a semiflexible chain 10, 22, 23, 24] and may 
result in the collapse of the PE chain . We have per- 
formed a calculation of the distribution function taking 
into account the attractive interaction due to counterion 
fluctuations and observed the signature of collapse. In 
general, one is able to observe the collapse of a charged 



13 



chain for any strong enough attractive interaction which 
is in function of the distance between monomers. We are 
currently investigating these effects [i^. 
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APPENDIX A: EXPANSION OF THE COULOMB 
INTERACTION 

In this Appendix, the expansion of the Coulomb self- 
interaction for a charged chain that is slightly deformed 
about the rodlike configuration is derived. Using Fourier 
representation of the screened Coulomb interaction we 
find 



[ ^ 

J (27r 



d^k 47r 



f fc2 _^ ^2 



ikz 



ds I ds' exp < zk^ • / du Si{u) I du a.{u)'^ -\- ikz{s' — s) 



1 f r' 

1 — - I kx • y du a.{u) 



ikz 
~2' 



du 8i{uf + O(a^) 



(Al) 



APPENDIX B: GENERAL FORM OF THE 
ENERGY ASSOCIATED WITH DISTORTIONS 
OF A ROD-LIKE POLYELECTROLYTE 
STIFFENED ENTIRELY BY UNSCREENED 
COULOMB INTERACTIONS 

Suppose that the tangent vector t{s) is as given by 
Eq. (@|). Suppose, also, that the interaction energy of 
the rod-like polyelectrolyte is given by 

£ = ^ [ ds j ds'V(v{s)-v{s')). (Bl) 

Making use of Eq. and the relationship 

v{s) = / t{si)dsi, (B2) 

and expanding to the resulting expression for the energy 
£ to second order in the deviation, a(s), from a straight 
line in th z-direction, we find that the energy is given by 



and 



£ = 
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and is the greater (smaller) of s and s' . If the 

charges on the polyelectrolyte interact via the unscreened 
Coulomb potential, then the interactions leading to the 
kernel C{s^s') are straightforward, and we find 
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Note that the kernel jC{s^s') is equal to zero whenever 
either one of the arguments is equal to or L. 

The kernel JC is neither local nor is it translationally 
invariant. However, if we are interested in what happens 
in the vicinity of 5 = 0, we can simplify jC and, as a 
consequence, JC. When both s and s' are much smaller 
than L, we can replace C as given by Eq. (|B6|) by 
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(B7) 



The kernel JC that results from this new C via Eq. (|B4|) is 
still non-local and is not translationally invariant. How- 
ever, it is possible to obtain, by inspection, its eigenfunc- 
tions and their associated eigenvalues. 

We assume that the minimum spacing between adja- 
cent charges, 6, is small, and we assume that we can 
retain those terms that are leading order in ratios of b to 



14 



other lengths in the system. Making use of these assump- 
tions, we are able to obtain, by inspection, results for the 
eigenfunctions and eigenvalues of the energy operator /C. 

We start by noting that the convolution of the energy 
kernel in Eq. (|B7|) with the function is equal to 

5^|21n^+p^ aP-Hn{l - a)da- 



P 



I aP-^ \n(a - l)da — ^ -\ 

Ji p^l p) 



(B8) 



The assumption underlying the calculations leading to 
(|B8|) is that the integrand is fundamentally convergent. 
That is, we ignore the possibility that the integrand 
yields non-integrable divergences anywhere in the region 
of integration. We now make use of the following rela- 
tions 



a' 



ln(l — a)da = —7 — ilj{p + 1), 
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where '^(x), the digamma function, is the logarithmic 
derivative of the gamma function: 
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This means that the "eigenvalue" associated with the 
eigenfunction is given by 



K{p) = 27 + V^(l + p) + V^(-p) + 



1 



p + 1 p 



1. (B12) 



Figure El display the function A(p) as a function of the 
real argument p. Note that the function is symmetric 




FIG. 17: the function A(p) as a function of real p. 

about p = 1/2. It can be readily shown that the function 
K{p) is real \i p = — l/2±z/c, with k real. The proper 
eigenfunctions and eigenvalues are associated with just 



such values of p. In fact, we can choose for eigenfunctions 
of the operator 



Us) 



s sin 



(fcln^ + 0(fc)). 



(B13) 



If we then require that the derivative of the eigenfunc- 
tion ik{s) is zero at the boundary (5 = 6), then (j){k) = 
arctan2/c. It can be shown that the integrations leading 
to the eigenvalue are all convergent, and the eigenvalue 
has the form A(— 1/2 + ik). A plot of this eigenvalue as 
a function of the parameter k is displayed in Figure ITSl 
This eigenvalue is an even function of k. As it turns out. 




FIG. 18: The function A(A:) = A(-l/2 + ik), for k > 0. 

an excellent numerical approximation to 

X{k) = A(-l/2 + i/c), (B14) 

is given by 

X{k) ^ 0.2274112 + In [l + 0.929A:^] (B15) 

The spacing of the eigenvalues is determined by the al- 
lowed values of the parameter k in (|B14|) . If we assume 
that the eigenfunctions pass through zero at some large 
value of s, or that they have a zero slope there, then 
it is straightforward to show that the lowest eigenvalues 
will be associated with equally-spaced /c's. Figure 01 con- 
tains a comparison between the eigenvalues given by Eq. 
(|B14|) with equally-spaced /c's and the numerically calcu- 
lated eigenvalues of the interaction operator with matrix 
elements in the cosine basis set as shown in Eq. 

As an additional check on the validity of the results 
presented here, we compare the eigenfunctions in Eq. 
(|B13|) to the eigenvectors of matrix in Eq. ©. The 
eigenvectors are plotted in real space. We find that for 
low eigenvalues, the expressions in Eq. (|B13|) provide an 
excellent match to the results of numerical calculations. 
Figure displays a comparison between an eigenfunction 
as given by Eq. ()B13|) and the numerically-determined 
eigenvector associated with the sixth-lowest eigenvalue of 
the discrete version of the kernel. The value of k in the 
formula for the eigenfunction was adjusted with the use 
of a least-squares procedure. The eigenfunctions have 
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zero slope at the boundary of the region associated with FIG. 19: The values of the parameter k. 

the smallest value of s. The quantity s does not range 
over the entire interval, from to 1, because the analytic 
eigenfunctions generated here are expected to be accu- 
rate only in the regime 5 <C where L has been set 
equal to one in the case at hand. 

The least-squares fitting of the analytical eigenfunc- 
tions to numerical eigenvectors leads to a set of values 
for the parameter, k. Figure El displays the values of 
/c, plotted against the index of the eigenvalues. As indi- 
cated by the straight line drawn through them, the /c's 
are approximately equally-spaced. 
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